Entanglement dynamics of two-qubit system in different types of noisy 

channels * 



C. J. Shant, J. B. Liu, W. W. Cheng, T. K. Liu*, Y. X. Huang, and H. Li 
College of Physics and Electronic Science, 
Hubei Normal University, Huangshi 435002, China 
(Dated: August 27, 2008) 

Abstract 

In this paper, we study entanglement dynamics of a two-qubit extended Werner-like state locally interact- 
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may exhibit entanglement sudden death (ESD) and/or be decreased for the critical probability at which the 
entanglement disappear. 

PACS numbers: 03.65.Ud, 03.67.Mn, 75.10.Pq 

Keywords: entanglement dynamics, entanglement sudden death, purity 



* Supported by the Natural Science Foundation of Hubei Province, China under Grant No 2006ABA055, and the 

Postgraduate Programme of Hubei Normal University under Grant No 2007D20. 
^ E-mail: scjll22@163.com 

Corresponding author. E-mail: tkliuhs@163.com 

1 



I. INTRODUCTION 



In the last two decades, entanglement attracts much attention due to the powerful applications 
in quantum information process and quantum computing In order to realize quantum infor- 
mation process, great effort has been devoted to studying and characterizing the entanglement 
in cavity QED and spin systemst^'^' schemes. Therefore, it is of increasingly importance to 
understand entanglement behaviors of quantum system in realistic situations, where the system un- 
avoidably loses its coherence due to interactions with the environment. In this context, a peculiar 
dynamical feature of entanglement was discovered that the entanglement can vanish completely 
in a finite time, in striking contrast with decoherence of its individual constituent that decays only 
asymptotically. Such a surprising phenomenon was termed entanglement sudden death (ESD)[^1. 
Because of its intrinsic and practical interests, ESD has attracted much attention in theory'^^^^i 
and confirmed experimentally '^'^l for the case of two qubits. The ESD phenomenon illustrates the 
fact that the global behavior of an entangled system may be markedly different from the individual 
and local behavior of its constituents. '^^1 

From a practical point of view, an entangled state undergoing ESD would put a limitation on 
the time of its application in practice since it is less robust than one without ESD. Hence, to find 
various conditions under which the ESD occurs seems to be very necessary. It has shown that ESD 
is sensitive to the type of initial entanglement, i.e., depending on the type of initial state, the en- 
tanglement may or may not exhibit ESD.l^'^l Besides the initial condition of an entangled state, the 
environment is another decided factor responsible for the dynamical behaviors of entanglement. 
In Ref. [14], the authors studied entanglement dynamics in three types of noisy channels. It was 
shown that for the same entangled states, their entanglement may exhibit completely different be- 
haviors involving the appearance of ESD in different types of environments. However, they only 
considered the pure case of initial entangled states. In this paper, we shall study the entanglement 
dynamics of an entangled state subject to various noisy channels by paying more addition to the 
initial condition of the state. We show that the purity of initial entangled state is closely related 
to the entanglement robustness in each noisy channel. That is, if the initial entangled state is pre- 
pared in mixed instead of pure form, the state may exhibit ESD and/or be decreased for the critical 
probability at which the entanglement disappear. 

To quantify the entanglement of a two-qubit system we adopt Wootters' concurrence'^^!. The 
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concurrence C{pab) for the density matrix pab of a two-qubit system AB is defined as 

C(Pab) = max{0, a/A^- a/A^- a/A^- a/aI}, (1) 



where Aj are the eigenvalues of the matrix C, = p{ciy ay)p*{ay ® ay) arranged in decreasing 
order. Here cr^*"^^ are the y-Pauli matrix acting on qubit A (B) and p* is the complex conjugation 
of p in the standard (computational) basis. For the separate state C = whereas C = 1 for 
maximally entangled state. 

For the initial state of qubit-pair AB, instead of Bell-like and Werner states '^^J, we shall con- 
sider the following extended Werner-like state 

Pab{0) = r \<^)abab ("^I + ^^^AB, (2) 
with r the purity of the initial state of qubits AB, Iab the 4x4 identity matrix and 

|$)^^ = (sine|00) + cos^|ll))^^, (3) 

the Bell-like state. Obviously, the state in Eq. ^ reduces to the standard Werner state when 
9 = Tx /A and to Bell-like pure state when r = 1. By dealing with the above extended Werner- like 
state, we are able to study the effect of mixedness of the initial entangled state. Both the Bell- 
like state and Werner state, and so the extended Werner-like state, belong to the so-called X-class 
state'^^1 whose density matrix is of the form 

( X ^ v\ 
y u 
M* 2 
\v* w / 

with X, ?/, 2;, w real positive and m, v complex quantities. The X-class states have the property that 
the corresponding two-qubit density matrix preserves the X-form during the system evolution. 
For the X-state dH), the concurrence can be derived as 

C{pab) = 2max{0, \u\ — a/xw, |f | — ^/yz}. (5) 

We consider three paradigmatic types of noisy channels, i.e., amplitude damping channel, phase 
damping channel and depolarization channel. Each one of the two qubits A and B individually 
coupled to its own noisy environment, implying no any interaction, direct or indirect, between 



Pab 



(4) 
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them. The dynamics of each qubit is governed by a master equation that gives rise to a completely 
positive trace-preserving map £i (with i = A, B) describing the evolution as pi = £ipi{0), where 
Pi{0) and Pi are, respectively, the initial and evolved reduced states of the i-th subsystem. In the 
following, we shall consider the entanglement dynamics of the state © in the three noisy channels, 
respectively. 

II. AMPLITUDE DAMPING CHANNEL 

The first noisy channel we consider is the amplitude damping (AD) channel which can be 
represented via the Kraus representation ast^'^"^! 

Sf^'p, = EopiEl + EipiEl (6) 

with Eq = |0) (0| + y/1 — p |1) (1| and Ei = ^/p\0) (1| being its Kraus operators, p = p(t) = 
1 — e~^T* is the probability of the qubit exchanging a quantum with the bath at time t, and 7 is the 
zero-temperature dissipation rate. After time-evolving, the initial state of qubits A, B in Eq. (|2]) 
evolves into pab (p) which has the X-form (g]) with the matrix elements are given as 

X = ^-^(1 + r(cos^(6')p^ + sin^(6')), 
1 — v 

y = z = — - — (1 — + r cos^{9)p{l — p), 

w = {^-^ + rcos^{e)){l - pf, 
V = r sinO cos6{l — p), 

u = 0. (7) 

By virtue of Eq. Q, we can get the concurrence C{pab{p)) of Pab{p) as 

C{pab{p)) = 2max{0, \v\ - (8) 

If the initial state of qubits A, B is pure, i.e., r = 1, from the relation \v\ — ^/yz = we can 
get the critical probability p^ at which the entanglement disappear as Pc = For 9 < tc/A, 

Pc is always smaller than 1, meaning that the entanglement disappears before the steady state is 
asymptotically reached. '^^^ Thus < 7r/4 is the condition for the occurrence of ESD when the 
initial state of A, B is pure, i.e. r = 1. The relations between C{pab{p)) and p, 9 are plotted 
in FIG. 1 for r = 1. However, we shall pay more attention to the case of r < 1, i.e., the initial 
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FIG. 1: Cab = C{pab{p)) as functions of 9 and p for the case of r = 1. The BSD appears for 6 < 
7r/4 + kTT. 

entangled state is prepared mixed form. In this case, the critical probability pc = ^^ir^'co^^g^^-t+i^ ' 
which is related to both the degree and purity of initial entanglement in terms of 9 and r. The 
condition of appearance of ESD is derived as | sin 9 cos 9 \ — cos^ 9 < ^ ( ^ — 1 ) , from which we can 
see that the range of 9 for appearing ESD depends on the value of r. As an example, for r = 0.7 
we plot in FIG. 2 the concurrence C{pab{p)) as functions of p and 9, where we can see that ESD 
occurs for all the possible values of 9, in striking contrast with the case of r = 1 in FIG. 1 . Hence, 
ESD phenomenon is closely related to the purity of initial entangled state. 



III. PHASE DAMPING CHANNEL 



In this section, we consider the phase damping channel (PD) in which there is loss of quantum 
information with probability p, but without any energy exchange. It is defined as^^^J 

£^P^ = (1 -P)P^ +P(|0)(0|p.|0)(0| + |l)(l|p.|l)(l|). (9) 

After time-evolving, the diagonal terms of initial state in Eq. (O of qubits A, B remains the same, 
whereas the off-diagonal ones are multiplied by (1 — pY- By virtue of Eq. ([5]), we can get the 
concurrence C{pab{p)) of Pab{p) as 

C{pab{p)) = 2max{0, \v\ - (10) 
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FIG. 2: Cab = C{pab{p)) as functions of 9 and p for the case of r = 0.7. The BSD appears for all the 
possible values of 9. 

with y = z = ^ and v = r sin6'cos6'(l — p^. Obviously, if the initial state of qubits AB is 
pure, i.e., r = 1, we have y = z = 0, thus C{pab{p)) = 2 sin6' cos6'(l — p)^ > with equality 
hold for p = 1 implying non-existence of ESD for all the possible values of 0. However, if the 
initial state of qubits AB is mixed the situation will be different. From the relation |t>| — ^Jyz = 
we can get the critical probability pc sls Pc = 1 — ^/i^^^^^^T^- For ^ < 1j Pc is always smaller 
than 1, meaning that in the PD channel the entangled state always suffers ESD if the state is 
prepared in the mixed form. To make a comparison, in FIG. 3 and FIG. 4, we plot C{pab{p)) as 
functions of 6 and p for r = 1 and r = 0.7, respectively. Hence, the close relation between ESD 
and the purity of initial entangled state can also be obtained in the PD channel. 



IV. DEPOLARIZING CHANNEL 



Next, we consider entanglement dynamics in the depolarizing (D) channel. The D channel 
represents the situation where the z-th qubit remains untouched with probability 1 — j9, or is 
depolarized-meaning that its state is taken to the maximally mixed state-with probability p. It 
can be expressed as'^^' 

£^p, = {l-p)Pi+p^, (11) 

with I is the identity operator. 
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FIG. 3: Cab = C{pab{p)) as functions of 9 and p for the case of r = 1. The BSD does not appear for all 
the possible values of 9. 




FIG. 4: Cab = C{pab{p)) as functions of 9 and p for the case of r = 0.7. The BSD instead appears for 
all the possible values of 9. 

After time-evolving, the density matrix Q evolves into Pab{p) which remains the X form dH) 
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FIG. 5: Cab = C{pab{p)) as functions of 9 and p for the case of r = 1. The BSD happens for all the 
possible values of 9. 

with the matrix elements are given as 

X = (l-|n^ + rsin^^)+Kl-|)i^ + (|f(i^ + rcos^^), 

y = ^ = ^(i-p+Y) + f(i-|)(i + r), 

w = (l-fH^ + rcos^^)+p(l-|)i^ + (|)^(i^ + rsin^^), 
V = r sm9 cos9{l — pY, 

u = 0. (12) 

By virtue of Eq. ([5]), we can get the concurrence C{pab{p)) of Pab{p) as 

C{pab{p)) = 2max{0, \v\ - (13) 

For r = 1 the state (|2l) always undergoes ESD for all the possible values of 6 in the time evolution 
process as shown in FIG. 5, where we plot C{pab{p)) as functions of p and 6. It is known that 
ESD is a representation for the fragility of entanglement. If we take the mixedness of initial 
entanglement into account, the robustness of entanglement will be reduced further. In FIG. 6 we 
plot C{pab{p)) as function of p for different r, where we can see that the critical probability will 
be decreased with r. 
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FIG. 6: C{pab{p)) as function of p for 6* = 7r/4 for various r. The critical probability decreased with r. 
V. CONCLUSIONS 

In conclusion, we have studied the entanglement evolution of a two-qubit entangled state sub- 
ject to independent environments, i.e., the amplitude damping, phasing damping and depolarizing 
channels. The initial entangled state is prepared in an extended Werner-like form thus we can 
investigate the relations between the purity of the entangled state and its dynamical behaviors. We 
find that the dynamical behaviors of the two-qubit entanglement are closely related to its purity 
as well as the types of noisy channels. The mixedness has direct impacts on the robustness of 
entanglement in the sense that it can result in ESD and/or decrease the critical probability. The 
understanding of entanglement behaviors in realistic system is a precondition for its application in 
practice, thus studying the conditions which may influence the entanglement dynamics in various 
situations prove very important and necessary. 
We thank Z-X Man for his reading of the manuscript. 
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